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Abstract 

Treating the quark and diquark as elementary particles, the Bethe-Salpeter equation 
for the nucleon is solved numerically The dependence of the mass on the diquark mass and 
on the coupling constants is investigated. The resulting relativistic quark-diquark-nucleon 
vertex is presented and discussed. 
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1 Introduction 

The quark-diquark picture has been successful in describing a variety of nucleon properties 
[jl|. It has been used to relate such different phenomena as the A-N mass difference 0], 
the neutron formfactor and the ratios of structure functions in deep inelastic scattering 
0, ^. In most applications one uses purely phenomenological quark-diquark wave functions. 
In order to actually calculate the dynamics of the nucleon, some authors have solved the Bethe- 
Salpeter equation for the quark-diquark system, with a local (quenched or static) quark-diquark 
interaction [Q. The integral equation is then separable and can be solved analytically. The 
exchanged quark does not propagate in this approximation. 

In this letter, the quark-diquark-nucleon (qdN) vertex is calculated, using the Bethe-Salpeter 
equation (BSE) for the quark-diquark system, with a propagating quark exchange potential. 
The properly normalized qdN vertex is shown to lead to the correct nucleon charges. Similar 
calculations, without inclusion of vector diquarks, have been presented in ^. 

The three-quark problem with a quark exchange potential has recently been solved by Ishii, 
Bentz and Yazaki [Q . However, they did not handle the problem of normalization and did not 
present results for the vertices (see our figure (H)), which are important e. g. for the calculation 
of formfactors and structure functions ^. 

2 The model 

We start from the following model Lagrangian: 

+ I g,^^C-^^\\x* - ig^MCrWx + g.^^ C-^l^r^^X^; - gl^^^Cr^^^x'^, 

where ■?/' is the quark field and rnq is the (constituent) quark mass. The scalar and vector 
diquark fields are denoted by x ^-iid x^ ■ The vector diquark field-strength tensor is denoted by 
F^^ = d^x^ ~(^^X^- The isospin index of the vector diquark is i = — , 0, +, and r~ = (1 — r^)/2, 
r° = l/\/2 and r+ = (1 + r^)/2 are the corresponding combinations of isospin SU{2) matrices 
for the quarks. The strengths of the quark couplings to the scalar and axial vector diquarks 
are gs and gy, respectively. The charge conjugation operator is C = i7^7°. 

The Bethe-Salpeter amplitude of the quark-diquark system has scalar and vector compo- 
nents (j){p, P) and (j)'^{p, P). Their mixing in the BSE is shown in fig. (0). The total and relative 
momenta in terms of the quark momentum pa and the diquark momentum pb are: 

P = VbPa - VaPb (2) 

P = Pa+ Pb, (3) 

where rja = nia/irna + rrib) and r^b = mb/{ma + rrib) with the quark mass nia = nig, and m^ 
is either one of the diquark masses rris or m„. We redefine the Bethe-Salpeter amplitudes as 
follows: 

4>{p) = [zS{pa)tS{pb)]-'Hp) (4) 

Up) = [^S{Pa)^S^^'{Pb)^'np), (5) 
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Figure 1: The mixing of scalar and vector 0^ quark-diquark-nucleon vertices in the Bethe Salpeter 
equation. 



thereby removing the external quark and diquark propagators: 
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The amphtudes (^^) satisfy the following equations 
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Here, (j)'^ is defined such that it involves the isospin " + " component of the diquark in the 
proton (the component with two up quarks in it). The vector-diquarks in the proton can have 
isospin " + " (up- up) or "0" (up-down) components. The corresponding vertices are related 
through isospin symmetry: 0^+ = —^/2(f)^ = (p^. The isospin " — " and "0" vertices in the 
neutron satisfy: 0"" = —^^cp^f* = —(p^. The momentum of the exchanged quark is denoted by 
Q- 



3 Structure of the solution 

The equations (|T^) describe an on-shell nucleon, with its spin |. Therefore the vertex functions 
(b and cA'^ can be written as: 



^f/(P, S) 



(13) 



0M = ^^f/(P,^), (14) 

where U{P, S) is the nucleon spinor and \1/ and \l/'^ are Dirac matrices. Multiplying the equations 
([T2|) by U{P, S) and summing over the nucleon spin orientation S, it is seen that only the 
projections: 

X = ^A+ (15) 

occur, where A"*" = (If + M)/{2M) is the positive energy projection operator. The on-shell 
vertex functions x ^i-iid x^ ^-^e Dirac operators that satisfy: 

X = XA+ (17) 

X'^ = X''A+- (18) 

They depend on the total 4-momentum P of the nucleon and the relative 4-momentum p of 
quark and diquark. Introducing the transverse relative momentum pi^ = p^ — P^-pf , one can 
make the following decomposition into scalar functions: 

for the scalar diquarks { p ■ a ■ P = p^a^uP'^ ), and: 

x' = v.jYt (i + I^) + ^^^'P' (^ "^ m) + ^^°^ 

M2 



^' (Me'^"^'^7,PpP, + M^^^'a'^'pT, + ^''p ■ a ■ PP^) + 
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V,^'P^[^T-j^p-a-P 

for the vector diquarks. The eight scalar functions 5*1,2 and V1...6, depend only on the scalars 
p ■ p, p ■ P and P ■ P = M^ (which is fixed). In the static limit, the vertices do not depend on 
the relative momentum. Therefore only the functions Si, V2 and V4 survive in this limit. In 
the non-relativistic limit, these functions should dominate, since the remaining functions are 
multiplied by positive powers of the relative 3-momentum. 

4 Bethe-Salpeter norm and nucleon charge 

From the inhomogeneous BSE, one can derive a normalization condition for the homogeneous 
BSE ||lO|. Defining the quark and diquark parts of the norm in the S(calar) and V(ector) 



channels: 



N,s ^ ^ I -0^,ltrxS{pa)^S{pa)S{p,)x (21) 

N,s = ' I 0^.lt''xSiPa)^^S\p,)x (22) 



N,v ^ ^ i ^\trrS{pa)^S{pa)S{p,Ux' 



N, 



dv = t I -0^^hrrSa,{Pb)j^{-2P-p,g^''+p^,P^ + P^pl)S,p{p,)x^, 



(23) 
(24) 



the BS-norm of the nucleoli can be expressed as a hnear combination of these one-loop integrals: 



r^:N,s + r]tN,s + -KN.v + vlN^v) = 1- 



(25) 



Calculating the electric formfactors of proton and neutron at zero momentum transfer, one 
finds for the charges of proton and neutron: 






P - ^Ngs + -Nds + -Ndv, 



l^r l.r ^.T ^.r 



(26) 

(27) 



If the vertex functions (|l^, ^OD are normalized according to eq. (|25|) , then these charges should 
be 1 and 0. 



5 Results and discussion 
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Figure 2: The mass M and binding energy E = M — nia — NsUis — N^m^ of the nucleon as a 
function of the scalar diquark mass, with m^ = ms-|-0.2 GeV and rua = rriq = 0.45 GeV . Dashed 
line: gs = 2.17, g^ = 1.64, in accordance with [|Tl], |I2|. Solid line: gg = 2.0, g^ = 1.51. 



In the rest frame of the nucleon, the scalar vertex functions depend on \p\ and pq. By taking 



appropriate traces of the BSE (|T^), one obtains coupled integral equations for the eight scalar 
vertex functions. 

The integrand of these kernels depend on \p\, pq and on the angle 6, defined hy p ■ jf = 
IpIIi^Icos^. This last dependence enters only through the potential and the corresponding 
integral can be done independently of the form of the eight scalar vertex functions. The 
integrands do not depend on the azimuthal angle 0. Poles in the propagators are avoided 
by performing a Wick- rotation: p° -^ ip^ . Now, one can start with eight arbitrarily chosen 
functions of 2 variables that are defined on some Nq ® N^ grid. Applying the BSE repeatedly 
one finds a stable solution if M has the correct value. In this sense an eigenvalue problem for 
the relativistic energy is solved. 



The parameters of the model are the scalar and vector diquark masses rris and ruy, the 
scalar and vector coupling constants gs and g^, the quark mass nia and the Euclidean cutoff A. 
We take m^ = 0.45 GeV for the quark mass, suggested by p and m^ — m^ = 0.2 GeV for the 
difference of the vector and scalar diquark masses, which reproduces the A — N mass difference 
[|]. The cutoff is taken to be A = 0.9 GeV. 

In figure (0) the mass M and binding energy E of the nucleon as a quark-diquark bound state 
are shown as a function of rUs- The binding energy is defined as £^ = M — rUa — Nsiris — N^m^, 
where A^^ = rj^Ngs + vt-^ds and N^ = ^{rj^Ngv + Vb^dv) are weighting factors for the scalar and 
vector channels. The nucleon mass is reached for a scalar diquark mass of about 0.6 GeV. This 
result however depends on the parameters of the model. In order to constrain these parameters 
further, it is necessary also to calculate other physical observables. Still it is interesting to 
see that the binding energy hardly depends on the constituent masses and is relatively low. 
Evidently in this picture most of the nucleon mass is already given by the sum of quark and 
diquark quasiparticle masses. 

If one puts the scalar coupling gg to zero by hand, it turns out that the potential in the vector 
channel is repulsive. However, if the scalar coupling is turned on, the potential is attractive in 
both the scalar and vector channels. 

The resulting Bethe-Salpeter vertex functions are plotted in figure (§) as functions of the 
absolute value of the relative 3-momentum p for different values of the relative energy po- This 
second variable does not have a non-relativistic analogue; the n.r. wave function is obtained 
from the relativistic vertex by integrating over the relative energy. The plotted functions include 
characteristic factors |j)|' where / = 0, 1, 2 correspond to s-, p- and d-waves. It is seen that the 
s-wave contributions Si, V\ and V4 dominate, but the p-wave contributions 5*2, V2, V5 and Ve 
are also important. The d-wave term V3 can be neglected. In the static limit, only the s-waves 
survive. 

The values of the norm in the various channels are: 

Ngs = 0.80, N^s = 0.76, A^.^. = 0.17, N^v = 0.14. (28) 

One observes that the conditions Qp = 1 and Qn = for the proton and neutron charge (see 
eqs. (p6|, |27|)) are very well realized. 

It is interesting to see that the charges carried by the S and V diquark channels in the 
proton, 0.79 and 0.20 respectively, differ substantially from the values obtained in the non- 
relativistic SU{2)spin ® SU{2)isospin limit where these charges are | for each channel. This 
result indicates that the scalar diquark component of the quark-diquark-nucleon vertex is more 
important than the vector diquark component. 

I would like to thank Wolfram Weise for fruitful discussions and careful reading of the 
manuscript. 
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Figure 3: The Bethe-Salpeter vertex functions as functions of the absolute value of the relative 
3-momentum p for three values of the quark-diquark relative energy p° = 0.05 GeV (solid), p^ = 
0.50 GeV (dashed) and p° = 1.61 GeV (dotted). 
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